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I. INTRODUCTION 

One of the most important frontiers in strong interaction physics is the study of the 
structure of the nucleon. Despite considerable experimental and theoretical progress made 
over the last forty years, there are still many unanswered questions. An example is the 
intensive debate which has continued over the spin structure of the nucleon, ever since the 
European Muon Collaboration (EMC) published their initial data on the spin structure 
function g x []]]. Traditionally, two types of observables related to the nucleon structure 
have been studied mostly extensively: elastic form factors and parton (quark and gluon) 
distributions. Electromagnetic form factors of the proton were first measured in the mid 
1950s @, and in recent years measurements of those of the neutron have been attempted and 
more are planned ||. Weak form factors are also being measured through parity- violating 
electron and neutrino scattering Q. On the other hand, the unpolarized quark and gluon 
distributions have been systematically probed in deep-inelastic scattering and Drell-Yan 
processes since the discovery of quarks at SLAC in the late 1960s. The polarized quark 
distributions have also been studied in a number of experiments in recent years ||, and 
more data on these are anticipated in the future from CERN, SLAC, HERA and RHIC. 

In this paper, we present a first detailed study of a new type of nucleon observable: the 
off-forward parton distribution (OFPD). The OFPDs generalize and interpolate between the 
ordinary parton distributions, measured for instance in deep-inelastic scattering, and the 
elastic form factors, and therefore contain rich structural information. There are no data 
so far on these distributions, and a relatively short theoretical history. To our knowledge, 
the OFPDs have come up independently in three different theoretical studies. In the late 
1980s, Geyer and collaborators || studied the relation between the Altarelli-Parisi evolution 
for parton distributions and the Brodsky-Lepage evolution for leading-twist meson wave 
functions. The "interpolating functions" introduced in Ref. [|| are essentially the OFPDs 
which we study in this paper. In the early 1990s, Jain and Ralston || studied hard processes 
involving hadron helicity flip, in terms of an "off-diagonal transition amplitude" involving 
off-forward matrix elements of bi-quark fields in the nucleon. It was shown that the integral 
of this amplitude over the quark four-momentum yielded elastic form factors. Recently, one 
of us H,[| introduced OFPDs in the study of the spin structure of the nucleon. The main 
observations in Refs. [^,[| were that the fractions of the spin carried by quarks and gluons can 
be determined from form factors of the QCD energy-momentum tensor, and that the latter 
can be extracted from the OFPDs. Furthermore, the deeply-virtual Compton scattering 
(DVCS) process was proposed as a practical way to measure the new distributions. 

From the point of view of parton physics in the infinite momentum frame, the OF- 
PDs have the following meaning: if a nucleon is moving with an infinite momentum in a 
particular direction, take out a parton with a certain fraction of the momentum, give it a 
four- momentum transfer A M , and insert it back into the nucleon. This is illustrated in Fig. 1. 
The OFPD is then the amplitude characterizing this process. On the other hand, from the 
point of view of elastic form factors, the moments of OFPDs are form factors of twist-two 
quark and gluon operators. For the spin-1 operators one has the ordinary electromagnetic 
and axial form factors, while for the spin-2 operators one has the form factors of the energy- 
momentum tensor. Because the form factors of the tensor contain information about the 
quark and gluon contributions to the nucleon angular momentum, the OFPDs can provide 



information on the fraction of the nucleon spin carried by quark orbital angular momentum 
— a subject of considerable current interest [|l(J . 

The OFPDs can be measured in diffractive processes in which the nucleon recoils elas- 
tically after receiving some momentum transfer. Moreover, one must have in the processes 
a hard, light-like momentum so that the parton light-cone correlations are selected. The 
simplest such process is DVCS [[J, in which a deeply- virtual photon, supplied by inelastic 
electron scattering, hits the nucleon and turns into a real, high energy photon. Such a 
process is easy to analyze theoretically and is similar to ordinary deep-inelastic scattering. 
More complicated processes include diffractive meson production, in which one must deal 
in addition with meson light-cone wave functions pT|-|I"5|j (these have recently been shown 
by Collins et al. |15| to be factorizable). The best experimental facility to carry out DVCS 
experiments is the proposed ELFE ||16|| . However, some studies can already be made at 
HERA (the HERMES collaboration), and at Jefferson Lab with a 6 GeV electron beam. 

In planning future DVCS experiments, it is important to have a theoretical estimate 
of the OFPDs. The purpose of the present study is to perform a first analysis of the 



OFPDs in the MIT bag model []T7|]. To be sure, the model has a number of well-known 
problems, including breaking of chiral symmetry and translational invariance, absence of 
explicit gluon degrees of freedom, etc. Nonetheless, it contains quarks; it predicts well the 
hadron spectrum; it gives reasonable initial input for quark distributions [p"8|-|20||; and it can 
describe the electromagnetic form factors fl9l , |2l|j of the nucleon. Other models of which we 



are aware have an equally long list of problems and do not seem to provide any obvious 
advantage for estimating the quark distributions. 

After firstly reviewing in Section II the definitions of the off-forward parton distributions 
and some of their general properties, we present in Section III the results of the bag model 
for their dependence on the various kinematic variables. The calculation accounts for the 
Lorentz boost and spectator quark effects. In Section IV we analyze the form factors of the 
energy-momentum tensor, and evaluate their t dependence. Finally, conclusions are noted 
in Section V, and possible extensions of this work outlined. 

II. BASICS OF OFF-FORWARD PARTON DISTRIBUTIONS 

In this section we review the definitions and model independent results for the OFPDs 
and their moments discussed in Ref. ||. Other definitions of OFPDs exist in the literature 
Pjnj, however, the definition introduced in Ref. M has a number of advantages, such as 
explicit hermiticity, and a simple connection with local operators and their form factors. We 
shall henceforth use the definition from Ref. ||. The OFPDs can be defined for both quarks 
and gluons, however, in this paper we focus primarily on quark distributions, since the 
description of the wave function of gluons in the nucleon is a much more difficult problem. 

To start, consider the bilocal operator ip(— An/2)£r M, 0(— \n/2), where A is a scalar 
parameter, ip is a quark field of a certain flavor, and T M = 7 M or 7^75. The light-like vector 
n^ is proportional to (1; 0, 0, —1), with a coefficient depending on the choice of coordinates. 
The gauge link C is along a straight line segment extending from one quark field to the 
other, which makes the bilocal operator gauge invariant. In the following, we work in the 
light-like gauge, A ■ n = 0, so that the gauge link can be ignored. 



One can now proceed to take the matrix element of the bilocal operator between the 
nucleon states of momenta P M and P' M = P M + A M , where A^ 1 is the four-momentum transfer. 
The matrix element must be expressible in terms of nucleon spinors, Dirac matrices, and the 
four-vectors P M , A M and n M . Since we are only interested in the leading-twist contributions 
which are proportioanl to P^ or P /At in the infinite momentum frame, we keep terms that 
are non- vanishing after multiplication by n^ [@||: 

'' e iXx {P'\^(-Xn/2)^(Xn/2)\P) = H(x,£,t) u{P')^u{P) 



2n 



+ E(x,U)u(P')^ W ^u(P) + --- , (la) 



! ' X e lXx (P'\i;(-\n/2)^ l5 i;(\n/2)\P) = H{x,£,t) u{P')^ l5 u{P) 



2vr 



+ E{x ) i ) t)u{P')^u{P) + ... , (lb) 



where t = A 2 and £ = — n ■ A, with u(P) the nucleon spinor, and the dots (■ • •) denote 
higher-twist contributions. It is possible to construct other Dirac structures that appear 
to be leading- twist, however, using Gordon identities and throwing away sub- leading terms 
one can always reduce these to the form in Eqs.(|]). The structures in Eqs.(|l]) are the 
same as those in the definition of the nucleon's elastic form factors. Examination of the 
helicity structure of quark-nucleon scattering shows that there are exactly four independent 
amplitudes. The chiral-even distributions H and H survive in the forward limit in which 
the nucleon helicity is conserved, while the chiral-odd distributions E and E arise from the 
nucleon helicity flip associated with a finite momentum transfer. 

The OFPDs are depicted graphically in Fig. 1, where k^ and fc' M are the four-momenta 
of the active partons. The physical meaning of the distributions becomes clearer if one 
introduces a conjugate light-like vector p^ of n^, with p-n = 1. Expanding P^ = (P + P') M /2 
and A M in terms of the vectors p^ and n M then gives: 

P M = jf + (M 2 /2K , (2a) 

A^ = -£ (p" - (M 2 /2K) + K, (2b) 

where M = P = M 2 — t/4, and the spatial components of P have been chosen along 
the z-direction. If we focus on the p M components of the momenta, the initial and final 
nucleons have longitudinal momenta (1 + ^/2)p^ and (1 — £/2)p M , and the outgoing and 
incoming quarks carry (x + C,/2)p fl and (x — £/2)p M , respectively. Since the nucleon cannot 
have negative longitudinal momentum, the limit on £ is obviously: 

< £ < 2 . (3) 

A more careful analysis using A^ > leads to the more stringent constraint 
< £ < \/—i/M. On the other hand, since quarks cannot carry more logitudinal 
momentum than the parent nucleon, one has the constraint on x: 

- 1 < x < 1 . (4) 



When x > £/2 both quark propagators in Fig. 1 represent quarks. When x < — £/2 they 
represent antiquarks. In these regions, the OFPDs are analogous to the usual parton dis- 
tributions. In the intermediate region, — £/2 < x < £/2, the quark propagators contain 
one quark and one antiquark, and here the distributions resemble a meson's wave function 
amplitude. 

The off-forward parton distributions display characteristics of both the forward parton 
distributions and nucleon form factors. In fact, in the limit of A^ 1 — * 0, one finds ||: 

H(x,0,0)=q(x), H(x,0,0) = Aq{x), (5) 

where q(x) and Aq(x) are the forward quark and quark helicity distributions, defined through 
similar light-cone correlations [[2l|. It must be pointed out that while the A^ — > limit is 
quite simple and natural for the OFPDs, it cannot be taken literally for the kinematics of the 
DVCS process, where a finite t-channel momentum transfer is essential to simultaneously 
maintain the initial photon deeply-virtual, and the final state photon real. 

On the other hand, the first moment of the off-forward distributions are related to the 
nucleon form factors by the following sum rules |7||||: 

l 
dxH(x,£,t) =F 1 (t) , (6a) 

i 

dxE(x,£,t) =F 2 (t) , (6b) 

i 

i __ 

dxH(x,£,t) =G A (t) , (6c) 

dxE(x,£,t) =G P {t) . (6d) 

i 

Here Fi(t) and F 2 (t) are the Dirac and Pauli form factors, and GU(£) and Gp(t) are the 
axial-vector and pseudo-scalar form factors of the nucleon. The x-integrated distributions 
on the left-hand-side of Eqs.(H) are m fact independent of £. The sum rules (||) provide 
important constraints on any model calculation of the OFPDs. 

Generalizing the sum rules (||), let us multiply Eg. ([la]) by x n ~ x and integrate x from —1 
to +1, which gives: 

n w n w . . . n ^(p'\$ iD»... iD ^-i 7 M»^|p) 
= £T B (£,t) u(P') i u{P) + E n {i,t) u(P') ) ^ V u{P), (7) 

where 

#»(£, t) = £ dx x n ~ l H{x, i, t), (8) 

and likewise for E n (£,t). The derivative D M is defined as 

D " = X - (p* - £") , (9) 



where 



i?" = 7?" + i^A", (10a) 

*D» = *d»- igA". (10b) 

The left-hand-side of Eq.(0) is an off-forward matrix element of the twist-two operator 

0$ v " M2 = i> i D {w ■ ■ • i D ^- 1 7 /x ™>-0 - traces, (11) 

where the braces {• ■ •} represent symmetrization of indices. On general grounds, a matrix 
element of O^ 1 "^ 2 is a sum of terms composed of form factors (which are functions of t only), 
appropriate Lorentz structures constructed from P , A M and the Dirac matrices, and nucleon 
spinors. The £ dependence in Eq. (|7|) arises only from contracting a vector A M in any of the 
Lorentz structures with the null vector n M . Therefore the £ dependence of the moments of 
the OFPDs is in the form of polynomials. To find the degree of the polynomials, notice 
that there are at most n contractions of n M with A M in the n-th moment. Thus H n (£, t) and 
E n (£,t) are polynomials of degree n in £, a result which is not obvious from the definition 
itself. The same considerations apply to the bilocal operator with 75 dependence. 

III. A BAG MODEL ESTIMATE OF OFF-FORWARD PARTON DISTRIBUTIONS 

In this section, we present a calculation of the OFPDs in a simple version of the MIT 
bag model |T7| . As mentioned in the Introduction, our choice of the MIT bag is based on 



the fact that it has quark degrees of freedom, and gives reasonable results for the electro- 
magnetic form factors |T9|j2T[] , as well as for polarized and unpolarized parton distributions 
|l8|-p0"|. However, being a model, it has a number of unwanted artifacts, such as a sharp 
boundary, absence of gluons, and the breaking of translational invariance and chiral sym- 
metry. Nonetheless, we believe that our results should provide a reasonable first guess of 
the unknown distributions at a low energy scale, 0(0.2 GeV 2 ). There are, of course, many 
other nucleon models on the market in which the OFPDs could be calculated, however, we 
see no clear reason that these models would be more reliable than that considered here. The 
issue of evolution of the distributions to higher energy scales will be addressed in a separate 
publication. 

When evaluating the OFPDs and form factors, it is convenient to work in the Breit 
frame, in which the initial and final momenta of the nucleon are: 

P, = (M; - At/2), P; = (M; £/2). (12) 

The t-channel momentum transfer squared becomes: 

t = - A* 2 = 4 (M 2 - M 2 ) . (13) 

Using Eqs.(2), we then have: 

p» = (1; 0, 0, 1)/(2M), n» = (1; 0, 0, -1)/M. (14) 

The variable £ in this frame is therefore related to the projection of A in the z direction, 

£ = -A z /M. (15) 
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For calculations in a model without exact translational invariance, the choice of frame is part 
of the model assumptions. In principle, a different result could be obtained if, for instance, 
one chose instead a frame where the initial nucleon was at rest. Nevertheless, we believe that 
the main features of our result will be weakly frame dependent, as in many other similar 
types of calculations. 

Recall that the coordinate space wave function of a quark in the rest frame of the MIT 
bag is given by: 

ij(r) = V^NR*(.ji eor ) ,) X , (16) 

\ia -r ji{e r) J A 

where R is the bag radius, eo = uj q /R is the quark energy, and u = 2.04 is lowest frequency 
solution of the bag eigenequation, tant^o = <^o/(l ~ ^o)- The functions j ,i are the spherical 
Bessel functions (r = \r\), and \ is the quark spinor. The normalization N is given by: 

N 2 = ^ (17) 

2iP (wb - 1) jg(wb) l ) 



The radius in the basic version of the bag model is given by the relation: RM = 4co>o Hl7| , |20 

Calculation of the OFPDs requires wave functions of a moving nucleon. One must there- 
fore boost the rest frame wave function ( |T6| ) to a frame moving with velocity v. Including 
the time dependence of the quark wave function explicitly, ip(t,r) = exp(— ieat) ip(r), the 
effect of a Lorentz boost on the wave function can be represented by [ pTJ| : 

i[)ff(t,r) = S (A.*?) if) ((£ — v • r) cosher; r + v- r (cosher — 1) — vt coshtu) 
= exp (—ie (t — v ■ r) cosh a;) 

xS'(A^) ip (r + v ■ r (coshtu — 1) — vt coshu;), (18) 

where 

b {Atf) = exp I — - — I = cosh — + v ■ a smh — . (19) 

\ Zi I £ Zi 

Here the rapidity u> is related to the velocity by uj = tanh -1 v, where v = \v\. In the Breit 
frame the velocity of the initial nucleon is v = — A /2M, so that 

= jj, smhw = — . (20) 

Here we basically treat the independent quarks in the bag as free particles, ignoring the fact 
that they are confined by the bag boundary, which again is a part of the model assumptions. 
In practical calculations, it will be more convenient to use a momentum space wave 
function, <f(k), which is simply related to the coordinate space wave function tp$(t,r) by a 
Fourier transformation: 

Mt, ?) = s ( A *) / 7^ ex P (-*(eo* - t ■ r)) <p(k), (21) 



where eo = (eo + k ■ v) cosily, k ± = k±, and k y = (eo^ + k\\) coshcj, with fe|| = {k ■ v) v. 
The momentum space wave function is given by: 

m = ^NX°^tl {k) )x, (22) 

where k = \k\, and the functions £0,1 are given by: 

. n\ JoM cos(kR) - j (kR) coscuo , n0 , 

to(k) = p— ; , (23a) 

Uq — k 2 R 2 

f (k s = 3o(kR) ji(u ) kR-j (uo) ji(kR) u 

One of the most important issues in calculating the off-forward matrix elements of single- 
particle operators in independent particle models like the bag is momentum conservation. In 
the following, we discuss this issue in some detail, which in the end will motivate the approach 
we take in performing the calculation. In independent particle models, strictly speaking the 
form factors of any one-body operator must be zero. Since the momentum transfer through 
the one-body operator affects the active quark only, the remaining spectator quarks maintain 
their original states. On the other hand, if a momentum transfer A is given to the whole 
nucleon, each of the quarks must receive a momentum transfer A/3, since before and after 
the interaction, the model nucleon must move as a whole. Thus, due to the momentum 
mismatch, the form factors must vanish. In the realistic situation, however, the nucleon 
wave function contains correlations. The momentum transfer injected to a single quark is in 
turn transferred through correlations to the other constituents, and asymptotically is equally 
shared among them. In independent particle models, however, these vital correlations needed 
for form factor calculations are missing. 

In the literature, several common approaches have been adopted to deal with this issue. 
In one approach, model wave functions with no definite center-of-mass momentum are used 
f21| , |T9|| . Form factors are calculated from a Fourier component of the single particle operator, 
0(A) = J d 3 f e tr ' A 0(f). In this type of calculation, the momentum transfer to the model 
nucleon and to the individual quarks is not well defined. If a boosted single-particle wave 
function is used, roughly speaking, the momentum transfer to each spectator quark is A/3. 
On the other hand, the momentum transfer through the active quark is approximately 
(1 — €q/M) A, which arises from the combined effects of boost and the action of the single- 
particle operator. (The relative sign of the two effects appears somewhat counterintuitive, 
but is nonetheless correct.) 

In principle, a better approach would be to use initial and final nucleon wave functions 
with definite momentum, which can be approximately obtained, for instance, through the 
Peierls-Yoccoz projection |22| or the center-of-mass freedom separation method [p3fl . In such 
calculations, all quarks share equally the momentum transfer to the nucleon. In particular, 
the active quark in which the single-particle operator acts is injected with a momentum 
A/3 only. Since the state of the two spectator quarks also changes from the initial to the 
final nucleon, the effective operator which induces such a transition is actually a three-body 
operator. 



A calculation of the OFPDs incorporating the effects of Lorentz boosts and projections 
is rather involved, and in practice not particularly illuminating. Instead, we consider a 
simpler alternative, by modifying the momentum transfer through the active quark in the 
approach of Ref. ||21|| . Namely, we let the effetive momentum transfer through the active 
quark be rjA, where rj is taken to be a parameter. As mentioned above, rj = 1 — €q/M in an 
unprojected calculation. Ultimately, rj in our calculation will be determined by fitting the 
electromagnetic form factors of the nucleon, but one can expect rj to be around 1/3, in the 
spirit of the momentum-projected calculation discussed above. 

The matrix element of the bilocal operator is calculated using the boosted wave function 
of the active quark, 



2M / — e lXx I d 3 r e iAr? ip_#{-\n/2 + r) V ^(Xn/2 + r) 
2M 



cosh a; J (2n) 3 



d ' k S (x - n_(k + + ^±)) [p\k') S(A_*) 7o r S(A,) <p(k)} 



(24) 



where T =%i or %i 75, and k' = \k'\, k' = k+ A. The effective momentum transfer A is given 
by: 



A = 7] 



A 



(25) 



coshu; 

Choosing for simplicity A y = 0, and using cylindrical coordinates to perform the k integra- 
tion, the (5-function in Eq.(f2~4|) reduces to a constraint on the ^-component of k: 

M 



A', 



x 



l-(coshcj-l)A2/t 

x = ( (2e n + A z ) cosh ijj + I A I sinh to 

2MV J ' ' 



2A X A, 



^(coshu; — 1' 



(26) 



Evaluating the expression in the braces in Eq. (|2~4f) explicitly for T =ti, and equating the spin 
independent components on both sides, leads to: 



H(x,Z,t) + -^E(x,£,t) 

Z 2 (t) UirN 2 R 6 



M 



dk± dip 
1 - (cosher - l)A 2 Jt J (2tt) 3 



x {t (k)t (k') 



+ 
+ 



k z cosh uj + 



k' z cosh uj 
2A 



2A 



t 



z k'-~& sinh 2 - 



fci 



to(fc) hjk') 
k' 



U) 



k ■ A sinh 2 
t 2 



h{k) to(fc') 



k ■ k' — A ■ ( k A 2 — A k z ) sinh uj 



k 
ti(k) hjk') 
kk' 



(27) 



where the effects of the spectator quarks are included in the factor Z(t) ||21|| : 

N 2 rR 



Z(t) 



cosh a; 



dr r z j e |A|r/M [j$(e r) + jf(e r) 



(28) 



If one compares the a y components in Eg. (|24j) , on the other hand, a different combination 
of H and E arises: 



H(x,£,t) + E{x,Z,t) 



2MM 



w V / l-(coshw-l)A2/t7 (2tt) 3 



f sinh cu . , . . , ,. 
x ( -^-toik) to(fc') 



I |A| 



A, 



2A 2 (A; x A) y . 2 w 
coshcj - sinh — 



A 



cosh uj 



2A Z (k x A) 
t A r 



sinh 



<jj 



t A x 
h(k) t (k') 



to(fc) ti(fc') 



A- 



+ 



(hA) s k-k'-2k 2 . 



A.., 



sinhtu 



*i(fc) ti(feQ 
kk' 



(29) 



Expressions for the individual functions if and £? are then obtained by solving Eqs.([27|) and 
([2|) directly. 

For the helicity dependent case V =7/1 75 again one can obtain two independent combina- 
tions of H and E by comparing different spin components in Eq.(^). Equating coefficients 
of a z leads to: 



2 1 



A? 



i?(x,£,t) 



A? 



-J(x,U) 



1MM 
2 1 / r dk ± dip 



x 



4M(M + M) 

Z\t) UttN 2 R 6 ) ' 2M — - / 

V / l-(cosha;-l)A2/ti (2vr) 3 

*o(*0 *o(fc' 



fci 



2A2 w 

cosh a; -\ smh — 

t 2 



k' H ^-^ sinh a; 

IAI 



to(fc) *l(fc') 

fc' 



+ 
+ 



A £* 

re, =; — smh 00 

IAI 



\2k z k' z — k ■ k'j cosh 
2A 7 



ti(fc) t (fc') 



re 



a; 



+ 



f 



(re, (£ + £') -A- A^ 2 )sinh : 



w 



kk' 



(30) 



while the a^ components give: 

A X A, 



2M(M + M) 
Z 2 (t) (4ttN 2 R [ 
' "2A T A 



2M 



smh — 



t 



dk± dip 
I - (cosher - l)A 2 Jt J (2vr) 3 

*o(fc) *o(*') + 



re. 



k' x %_E gjn^ w 



|A| 



to(fc) ti(fcQ 
re' 



+ 



fea: H =i — smh uj 

IAI 



ti(fc) t (fc') 
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(k x k' z + k z k' x ) cosh a; 



f 



k x (k + &)■£- A X P) sinh^l tl( ^ (fc/) | . (31) 



By fixing the bag radius to be R = Auq/M [[17],|(J, the model then has essentially one 



parameter, r\. This can be constrained by comparing the model predictions for the Sachs 
electric and magnetic form factors of the proton, 

G E (t)=F 1 (t) + ^F 2 (t), (32a) 

G M (t)=F 1 (t)+F 2 (t), (32b) 

with the available data |24|,|25|. The form factors can be obtained from Eqs.(^) and (|6b]) 



by integrating the H and E distributions in Eqs.(^Tj) and (^) directly. Note that the 
results are in fact equivalent to those given for Ge and Gm in coordinate space in Ref. 
| 2"T| . Furthermore, in a calculation for the proton with an SU(6) symmetric wave function, 



the combinations of the electric, baryonic, and axial charge factors appear in such a way 
as to reduce the final result to that for a single quark. In Fig. 2 we show the predicted t 
dependence of Ge for two values of the parameter 77 (77 = 0.35 and 0.55). The t dependence 
of the magnetic form factor Gm is shown in Fig. 3. In both cases the bag model results are 
in quite good agreement with the data. The small- 1£| data (including the charge radius and 
magnetic moment) do favor the larger value 77 = 0.55, while a better fit can be achieved 
at large \t\ with 77 = 0.35. Note that while G E (0) is independent of R and r/, the value of 
the magnetic moment G^(0) depends on the bag radius and rj. This shows that different 
prescriptions of calculating form factors within the model can give different answers even 
at small momentum transfers. This is just one of the artifacts of the explicit breaking of 
translational invariance in the bag model. 

Having fixed the model parameters, one can now calculate the individual OFPDs as a 
function of x and £, for different values of t. Again, assuming the SU(6) wave function for 
the proton, one multiplies the right-hand-side of Eq. ( |2T| ) by a factor 2 (1) and that of Eqs. 
(H|), fl30D and ( |31~D by a factor 4/3 (-1/3) to solve for the up (down) quark distributions. 
In Figs. 4(a) and (b) we first show the distributions at t — (and £ = 0) for — 1 < x < 1, 
for both the u and d quark flavors. Because the small- \t\ form factors are better described 
with a larger 77 value, for consistency we use here 77 = 0.55. Note that the distributions 
H and H in Fig. 4(a) for u and d quarks are just the forward unpolarized u(x) and d(x), 
and polarized Au(x) and Ad(x) parton distributions, respectively. Both H and H are in 
fact independent of 77. Furthermore, the first moments of H and H at t — explicitly 
satisfy the normalization conditions (equal to 2 and 1 for unpolarized u and d quarks, and 
4/3 gA and —1/3 5% where qa = 0.65 for a quark, in the polarized case, respectively). The 
tensor and pseudoscalar distributions E and E, however, shown in Fig. 4(b), do depend 
rather strongly on 77. The peak value of E for the u quark for instance would be ~ 6 for 
rj = 1 — €q/M = 0.75 (with R = 4uq/M). The first moment of E at t = 0, summed over 
both the u and d flavors, is equal to Gm(0) — 1 (as in Fig. 3), while the first moment of the 
flavor singlet E is Gp(0) ~ 6 for rj = 0.55. Note also that the calculated x distributions do 
not vanish entirely at x — 1, which simply reflects the fact that the initial and final nucleons 
are not good momentum eigenstates. However, the effect is only slightly noticeable for E, 
and is negligible for the other distributions. 
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In Figs. 5-10 the distributions for both the u and d quarks are shown as a function of 
x and £, for t = —1 and t = —2 GeV 2 . Throughout we use the value rj = 0.35 on account 
of the better agreement with the form factor data at large \t\ (see Figs. 2-3). One can see 
that the dependence on £ is quite weak. According to the discussion in the last section, 
this means that the form factors of the twist-two operators associated with the structure 
p 1 ...p n ~ 1 u{P') n f IXn u(P) dominate over other form factors. We do not have a simple 
explanation for this. The t dependence of OFPDs, however, is rather strong, as expected 
from the above form factor behavior. Experimentally, the most interesting region for DVCS 
is 1 < — t < 2 GeV 2 . For too small \t\, QED radiative effects mask the processes sensitive 
to the OFPDs. On the other hand, for too large \t\, the distributions become too small to 
be measurable. 

At asymptotically large \t\, the OFPDs can be analyzed using perturbative QCD. Again, 
it is the leading-twist light-cone wave function which determines the behavior of the helicity 
conserving distributions. We expect therefore H(x,£,t) and H(x,£,i) to fall like 1/t 2 as 
—t — ► oo. 

IV. FORM FACTORS OF THE ENERGY-MOMENTUM TENSOR 

In this section we review the role of the form factors of the QCD energy-momentum 
tensor played in the spin structure of the nucleon, and their relation with the OFPDs. In 
particular, we present the first model calculation of the form factors of the quark part of the 
energy-momentum tensor. 

Since the publication of the EMC measurement [|]] of the fraction of the proton's spin 
carried by quarks, the spin structure of the nucleon has been studied extensively in the 
literature. A deeper understanding of the problem, however, requires one to examine more 
closely the angular momentum operator in QCD. This can be written as a sum of quark and 
gluon contributions ||: 

J qcd = J q + J g , (33) 

where 



~jt q = J d 3 rfx T^ 



d 3 r 



-ijj^ip + ^ r x (-il5)ip 



2 



(34a) 



~f g = j d 3 r fx x if) . (34b) 

where summations over flavor and color indices are implicit, and T q and E x B are the 
quark and gluon momentum densities, respectively. The Dirac spin-matrix is denoted by U, 
and JJ is the covariant derivative. By analogy with the magnetic moment, the separate quark 
and gluon contributions to the nucleon spin can be obtained from the form factors of the 
momentum density, or equivalently the QCD energy-momentum tensor at zero momentum 
transfer. 

The symmetric, conserved, gauge-invariant energy-momentum tensor T^ u of QCD can 
be separated into quark and gluon components: 
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where the quark part is: 

Tf = - (iJj^W^ip + ^ 7 WS^v) , (36) 

and the gluon part is: 

Tf = ^F 2 - F^F^ . (37) 

Using Lorentz covariance and invariance under the discrete symmetries, one can expand the 
matrix elements of Tjj*i in terms of four form factors: 

{P , \T^\P) = u{P , )[A q , g {t)^T } + B q!g (t)P { "ia u ^A a /2M 

+ C q , 9 (t) (A" A" - g^t)/M + C q , g (t) g^M]u{P) , (38) 

where the braces {■ • •} on the superscripts denote symmetrization. Substituting the above 
into the nucleon matrix element of J q , g , one finds fractions of the nucleon spin carried by 
quarks, J q , and gluons, J g ||, 

J q *=l(\g(P) + B q , g (0J) , (39a) 

J q + J g = \, (39b) 

where in the rest frame of the nucleon 

J q J=(P\7 q jP), (40) 

with S the polarization vector of the nucleon. This role of the form factors of the total 



energy-momentum tensor has been first noted by Jaffe and Manohar |26 . 

Measuring the form factors of the energy- momentum tensor in practice is very difficult, 
partly because there is no fundamental probe which couples to them. (The graviton does, but 
only to the total tensor.) Because of asymptotic freedom, the form factors can be measured 
through deep-inelastic sum rules, as explained in Ref. |J. According to our definition of 
OFPDs, it is simple to show that: 

J dx xH(x, f , t) = A(t) + £ 2 C(t) , (41) 

J dx xE(x, f , t) = B(t) - £ 2 C(t) . (42) 

Combining these one obtains: 

f dxx(H(x,£,t) + E(x 1 £,tj)=A(t) + B(t) , (43) 

where the £ dependence drops out! 
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For the spin structure of the nucleon only small values of t are relevant. However, 
reaching small t in real experiments will be difficult; so some knowledge of the t dependence 
of the form factors is essential. In the remainder of this section, we will present a detailed 
calculation of the form factors of the (quark flavor-singlet part of the) energy-momentum 
tensor in the bag model. In principle, one can already obtain the form factors from the 
sum rules (|4lD-(|43|). However, here we will present an independent derivation, from which 
not only can one check the calculations on OFPDs, but also obtain the form factors not 
accessible from the above sum rules. 

The bag energy-momentum tensor is a sum of quark and empty-bag contributions. The 
quark part can be written: 



rpl-lU 
1 



- (ijjj {tM id^tfj + ipj^id^if) 



(44) 



where ip is the quark field in the bag. Using the bag's equation of motion, i $ i\) = S(r — R), 
it is easy to see that T^ B is traceless. From the general theorem proved in Ref. f27H , the 
bag quarks contribute 3/4 of the nucleon's mass. However, as we shall see below, the quarks 
in the bag carry all of the momentum of the nucleon. 

The form factors A,B,C and C are calculated from the matrix elements of T^ V B by 
taking various components of T^ B and comparing with the form in Eq. (|38|) . Taking A to 
be in the ^-direction, one finds from the T® l B components: 



+ 



M 



eo 



r\ t 



iai \M m l 



t (k)t (k') 



lt z K, 



kk 



'fhifyhik') sinhw 



+ (^t (k)ti(k') - jti(k)to(k')j cosh 



LO 



(45) 



In the case of no Lorentz boost, u — 0, and the combination A + B is unity at t — 0. 
According to Eq.( |39a| ), this represents the fact that the total angular momentum (spin plus 
orbital) of quarks in the bag is 1/2. In Fig. 11 this is illustrated by the dotted curve. With 
a Lorentz boost and a momentum transfer fraction to the quark rj < 1 — €q/M, the value of 
A(0) + -B(O) is no longer unity. This can be understood from the fact that the boosted bag 
wave function does not have the correct Lorentz symmetry. For a smaller value rj = 0.35 
(solid) one finds that A + B at t = is now w 0.5. However, as discussed in the previous 
section, the small-t form factors are more accurately described with a larger t], rj = 0.55 
(dashed) which gives A(0) + B(0) w 0.7. 

In dispersion theory, the t dependence of form factors is controlled by the nearest sin- 
gularities in t-channel. For A(t) + B(t), the quantum number of the channel is the exotic 
state with J PC = 1~ + . There is no conclusive evidence at the present time for the existence 
of resonances in this channel, although theoretical investigations indicate that hybrid qqg 
mesons could exist in the mass range of 1.3 to 1.9 GeV [23. If the dispersive behavior of 



A{t) + B{t) is dominated by large- mass resonances, then A + B will vary slowly with t, at 
least noticeably slower than the electromagnetic form factor. Although the bag calculation 
does indicate such a trend, the evidence is not strong. We suspect therefore that either the 
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t dependence of the bag prediction is not reliable, or multi-pion cuts in the form factor are 
important. Further study in this direction is called for. 

To obtain the individual form factors A, B, C, C, one can take other components of the 
tensor T^ Bag . The T®° Bag component gives: 

= 3 Z 2 (() (ircN'R") (^ - ±^\ j ^ {t«(k)t (k') + k-k> tl (k) tl (k')} , (46) 

where the factor of 3 is the number of valence quarks, and from the T® B (or T gB ) 
components one gets: 



t {k)t x {k') t x {k)t Q {k') 



coshw + v ' / ; [A|sinha;^. (47) 



k> k 

Finally, the T** Bag components give a fourth equation: 

C(t) = -3 Z\t) (47riV 2 i? 6 ) ^j-^- 3 (k Z + K)[^to(k)h(k') + ^hikMk')^. (48) 



Solving Eqs.(^5|)-(|48|), we plot in Fig. 12 the resulting form factors as a function of t for 
f] = 0.35. All form factors fall off monotonically as —t increases, with A(t) being the 
largest (and positive), and others relatively smaller (and negative). Note that for a radius 
R = 4u>o/M, the form factor A(0) = 1, which according to the definition of the form factors 
simply reflects the fact that all the momentum of the nucleon is carried by quarks. 

The trace condition on the quark part of the bag energy-momentum tensor implies that 
not all of the form factors are independent. In fact, it is easy to show that, 

t 3/ _ 

A (t) + m - 2 B(t)-—C(t) + AC(t) = 0. (49) 

We have checked that this relation is in fact satisfied to a very good approximation in the 
bag model, despite the violation of Lorentz symmetry after boosts. The above equation 
indicates that the normalization of the non- conserving form factor is fixed to C(0) = 1/4. 
It must be pointed out, however, that the above relation is not true in QCD because the 
renormalized quark part of the energy- momentum tensor has a trace anomaly p9| , |30" l. 



V. CONCLUSION 

In this paper we have presented a detailed study of a new class of nucleon observables, 
the off-forward parton distributions. The physical significance of the distriubtions has been 
explained in a partonic language, and their relations to form factors of twist-two operators 
made explicit. The £ dependence of the moments of the OFPDs is found to have a rather 
simple polynomial form. We have made the first model calculation of the distributions 
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using the MIT bag model, with the specific details of the model fixed by requiring one 
to reproduce the electromagnetic form factors of the proton, and the gross features of the 
parton distributions at a low energy scale. In relation to the spin structure of the nucleon, 
we have also studied the form factors of the energy-momentum tensor in the bag model, 
focusing on the t dependence in the range of to -2 GeV 2 . 

The model calculation includes the effects of Lorentz boosts of the quark wave function, 
and spectator quarks. The £ dependence of the distributions turns out to be extremely 
weak, indicating that the form factors of the twist-two operators associated with the Lorentz 
structure P Ml ■ ■ ■ P^"~ 1 u(P')'j fin u(P) are dominant. The t dependence of the OFPDs, on the 
other hand, is much stronger, and roughly follows that of the elastic form factors. Our results 
for the combination A + B of the tensor form factors provide the first concrete indication of 
their possible t dependence, which will be important for the extraction of the value at t = 0, 
since this gives the total angular momentum carried by quarks in the nucleon. 

An important issue to be addressed in future is the Q 2 evolution of the OFPDs. It is 
well known that the bag model calculations of the forward parton distributions are valid 
only at very low energy scales, 0(0.2) GeV 2 , therefore the OFPDs calculated in this paper 
cannot be used directly to estimate the experimental cross sections. The evolution equations 
for the off-forward distributions, which interpolate the Brodsky-Lepage |HJ and DGLAP 
32] evolution equations, have been derived in Refs. pj,!!,!^]. The technology for solving 



these new evolution equations is currently being developed. Apart from the evolution, the 
model calculations of the OFPDs could also be refined in future by including Lorentz boosts 
together with momentum projections. 

The most relevant experiment from which the off-forward distributions and tensor form 
factors can be extracted is deeply- virtual Compton scattering |||§. In practice, the DVCS 
process will be overwhelmed at small t by the Bethe-Heitler process, whose cross section 
has a QED infrared divergence at t — 0, so that one cannot go to too small t. One may get 
around this to some extent by subtracting the (calculable) Bethe-Heitler contribution, or by 



isolating the virtual Compton and Bethe-Heitler interference term |33| through single-spin 
asymmetry or combined lepton-antilepton scattering ||. Other processes which can also 
provide information on the OFPDs include diffract ive p and J/ip production |12-l5|1, which 



are sensitive to the off-forward gluon distribution. A series of dedicated experiments at 
suitable machines, such as HERMES, ELFE, or at Jefferson Lab, can in future explore these 
new distributions, and thus provide valuable information on the distribution of spin in the 
nucleon. 
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FIGURES 




FIG. 1. Off- forward parton distribution. P and P' are the initial and final state nucleon 
momenta, and k and k' are the active quark momenta. 
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FIG. 2. t dependence of the proton electric form factor Ge(£) for rj = 0.35 (solid) and rj = 0.55 
(dashed). The bag radius in this and subsequent figures is fixed by RM = 4ujq. The data are from 
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FIG. 3. t dependence of the proton magnetic form factor Gm(^) for rj 
rj = 0.55 (dashed). The data are from Ref. [pU2l 
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FIG. 4. Off-forward parton distributions at t = for the u and d flavors: (a) i7 and H, 
which correspond to the usual spin averaged and spin dependent distributions q{x) and Aq(x), 
respectively; (b) E and E for rj = 0.55. Note that E/2 is plotted. 
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FIG. 5. Off-forward parton distribution H (x, £, t) for the u quark, as a function of x and £, for 
(a) t = -1 GeV 2 and (b) t = -2 GeV 2 . 
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FIG. 6. As for Fig. 5, but for the d quark. 
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FIG. 7. Off-forward parton distribution H (x, £, t) for the u quark, for (a) t 
t = — 2 GeV 2 . The d quark distribution is obtained by multiplying by -1/4. 



-1 GeV 2 and (b) 



24 



0.6 



£ = -1 GeV 



2 




0.6 



Q.4 






0.2 



(b) t=-2GeV 
u quark 



2 




FIG. 8. Off-forward parton distribution E(x, £, t) for the w quark, for (a) t 
t = -2 GeV 2 . 



-1 GeV 2 and (b) 



25 



0.8 




0.8 



0.6 



H 0.4 

i 

0.2 



(b) t=-2Ge"~ 
d quark 



2 




FIG. 9. As for Fig. 8, but for the d quark. Note that it is — E(x,£,t) which is plotted. 



26 



1 .2 



£ = -1 GeV 



2 




7^ 



1 .2 



O.S 



0.4 



(b) t=-2GeV 
u quark 



2 




FIG. 10. Off-forward parton distribution E(x, £, t) for the u quark, for (a) t 
(b) t = —2 GeV 2 . The d quark distribution is obtained by multiplying by -1/4. 



-1 GeV 2 and 



27 



+ 








u. 



H 



o 



-t 



2 



FIG. 11. t dependence of the form factor A + B of the nucleon energy-momentum tensor, for 
the unboosted calculation with rj = 1 — eo/M (dotted), and for the boosted result with r/ = 0.55 
(dashed) and rj = 0.35 (solid). 
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FIG. 12. t dependence of the form factors A, B, C and C of the energy-momentum tensor, for 
rj = 0.35. 
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